Heavy Quarks and Leptons, Melbourne, 2008 



1 



Yus determination from kaon decays 

P. MassarottO 

Universita' degli studi di Napoli "Federico 11" and INFN Napoli, Naples, Italy 

This review presents the analysis of leptonic and semileptonic kaon decays data done by the 
FlaviaNet Kaon Working group, as described in Data include all recent results by BNL-E865, 
KLOE, KTeV, ISTRA+, and NA48. Experimental results are critically reviewed and combined, 
taking into account theoretical (both analytical and numerical) constraints on the semileptonic 
kaon form factors. We report on a very accurate determination of Vus as well as on many other 
tests of the SM which can be performed with leptonic and semileptonic kaon decays. 
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INTRODUCTION 



In the Standard Model, SM, transition rates of 
semileptonic processes such as — >■ uHv, with (u^ ) 
being a generic down (up) quark, can be computed 
with high accuracy in terms of the Fermi coupling 
Gf and the elements Vji of the Cabibbo-Kobayashi 
Maskawa (CKM) matrix. Measurements of the tran- 
sition rates provide therefore precise determinations 
of the fundamental SM couplings 
A detailed analysis of semileptonic decays offers also 
the possibility to set stringent constraints on new 
physics scenarios. While within the SM all (f — )■ uHi/ 
transitions are ruled by the same CKM coupling Vji 
(satisfying the unitarity condition J^k l^ifeP — 1) ^^'^ 
Gf is the same coupling appearing in the muon de- 
cay, this is not necessarily true beyond the SM. Setting 
bounds on the violations of CKM unitarity, violations 
of lepton universality, and deviations from the V — A 
structure, allows us to put significant constraints on 
various new-physics scenarios (or eventually find ev- 
idences of new physics). In the case of leptonic and 
semileptonic kaon decays these tests are particularly 
significant given the large amount of data recently 
collected by several experiments: BNL-E865, KLOE, 
KTeV, ISTRA+, and NA48. The analysis of these 
data provides precise determination of fundamental 
SM couplings, sets stringent SM test almost free from 
hadronic uncertainties, and finally can discriminate 
between new physics scenarios. The high statistical 
precision of measurements and the detailed informa- 
tion on kinematical distributions have provided sub- 
stantial progress on the theory side, in particular the 
theoretical error on hadronic form factors has been re- 
duced at the 1% level. 

The paper is organized as follows. First in Sec. |lT] 
we present fits to world data on the leading branch- 
ing ratios and lifetimes, for Kl, Ks, and mesons. 
Sec. mil summarizes the status of the knowledge of 
form factor slopes from K^^ decays. The physics re- 



TABLE I: Results of fit to Kl BRs and lifetime. S is the 
scale factor applied to the error in order to obtain chi^ — 
ndf. 



Parameter 


Value 


S 




0.4056(7) 


1.1 




0.2705(7) 


1.1 


BR(37r°) 


0.1951(9) 


1.2 


BR(7r+7r-7r°) 


0.1254(6) 


1.1 


BR(7r+7r-) 


1.997(7) xlO"^ 


1.1 


BR(27r°) 


8.64(4) xlO"* 


1.3 


BR(77) 


5.47(4) xlO"* 


1.1 


TL 


51.17(20) ns 


1.1 



*WWW access at www.lnf.infn.it/wg/vus/; for sake of com- 
pleteness and brevity for all references we refer to the Note 
written by the FlaviaNet Kaon WG [J • 



suits obtained are described in Sec. IIVI in particular 
the measurement of \Vusf+{0)\- 



II. EXPERIMENTAL DATA: BRS AND 
LIFETIME 

Numerous measurements of the principal kaon BRs, 
or of various ratios of these BRs, have been published 
recently. For the purposes of evaluating \Vusf+{0)\, 
these data can be used in a PDG-like fit to the BRs 
and lifetime. A detailed description to the fit proce- 
dure and the references of all experimental input used 
can be found in Ref. [H. 

For Kl the results are given in table HI while table llll 
gives the results for K^. 

For the Ks, the fit is dominated by the 
KLOE measurements of BR{Ks — irev) and of 
BR{tt+'k-)/BR{t:°'k°). These, together with the con- 
straint that the Ks BRs must add to unity, and 
the assumption of universal lepton couplings, com- 
pletely determine the Ks leading BRs. In particu- 
lar, BR(Xs ^ Tiev) = 7.046(91) x 10"''. For tks 
we use 8.958 xlO~-^^ s, where this is the non-CPT 
constrained fit value from the PDG, 

The fit takes into account the correlation between 
these values, as well as their dependence on the K^ 
lifetime. The world average value for t± is nom- 
inally quite precise; the 2006 PDG quotes t± = 
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TABLE II: Results of fit to BRs and lifetime. 



Parameter 


Value 


S 


BR(A-^2) 


63.57(11)% 


1.1 


BR(7r7r'') 


20.64(8)% 


1.1 


BR(7r7r7r) 


5.595(31)% 


1.0 


BR(A-e3) 


5.078(26)% 


1.2 


BR(7^p3) 


3.365(27)% 


1.7 


BR(7r7rV) 


1.750(26)% 


1.1 




12.384(19) ns 


1.7 



12.385(25) ns. However, the error is scaled by 2.1; 
the confidence level for the average is 0.17%. The two 
new measurements from KLOE agree with the 
PDG average, and give a smaller scale factor to the 
T± value. 



III. EXPERIMENTAL DATA: FORM 
FACTORS 



The hadronic K n matrix element of the vec- 
tor current is described by two form factors (FFs), 
f+{t) and /o(0, defined by {tt- (k) \sru\K° {p)) = 
{p + krf+{t) + {P ~ krf-{t) and f_{t) = 

'■^^^ (Mt) f+it)) where t = {p - kf . 
By construction, /o(0) = /+(0). In order to compute 
the phase space integrals we need experimental or the- 
oretical inputs to determinate the ^-dependence of FF. 
In principle, Chiral Perturbation Theory (ChPT) and 
Lattice QCD are useful tools to set theoretical con- 
straints. However, in practice the t-dependence of the 
FFs at present is better determined by measurements 
and by combining measurements and dispersion rela- 
tions. Many approaches have been used, and all have 
been described in detail in [ij. For K^z decays, re- 
cent measurements of the quadratic slope parameters 
of the vector form factor (A^, A" ) are available from 
KTeV, KLOE, ISTRA+, and NA48. 
The same collaborations recentely measured also the 
slope parameters (A^, A" , Aq) for K^^ decays Here we 
list only the averages of quadratic fit results for iCes 
and Kfj_j, slopes piip used to determine |T4s|/+(0). 
It is important to stress that the significance of the 
quadratic term in the vector form factor is strong for 
both Kf.^ (4.2(t) and i^^3 (3.6cr) fit to all data. 



TABLE III: Averages of quadratic fit results for K^z and 
slopes. 





Kl and K 


Measurements 


16 


xVndf 


54/13 (7 X IQ-'^) 


\'+ X 10^ 


24.9 ± 1.1 {S = 1.4) 


A^ X 10^ 


1.6 ±0.5 (S = 1.3) 


Ao X 10^ 


13.4 ±1.2 (5*= 1.9) 


P(AV,AT) 


-0.94 


p(A'+,Ao) 


+0.33 


p(A'|,Ao) 


-0.44 




0.15457(29) 




0.15892(30) 




0.10212(31) 




0.10507(32) 




+0.63 



mode 


|K.|/+(0) 


% err 


BR 


r 


A 


Int 


Kl -) 


■wev 


0.2163(6) 


0.28 


0.09 


0.19 


0.15 


0.09 


Kl -) 


TVflU 


0.2168(7) 


0.31 


0.10 


0.18 


0.15 


0.15 


Ks ^ 




0.2154(13) 


0.67 


0.65 


0.03 


0.15 


0.09 


K^ - 


> nev 


0.2173(8) 


0.39 


0.26 


0.09 


0.26 


0.09 


K^ - 




0.2176(11) 


0.51 


0.40 


0.09 


0.26 


0.15 


averaf 




0.2166(5) 













TABLE IV: Summary of jKis|/+(0) determination from 
all channels. 



IV. PHYSICS RESULTS 
A. Determination of |Kisl/+(0) and 

\Vus\/\V^cl\fK/U 

The value of |Kis|/+(0) has been determined from the 
decay rate of kaon semileptonic decays: 

r(A'£3(7)) = tI^^^^™ IK.|V+(o)' X (1) 

4(A+,o) [l + 5lui2)+€Sf 

using the world average values reported in previous 
sections for lifetimes, branching ratios and phase space 
integrals and the radiative and SU (2) breaking cor- 
rections discussed in [l[. The results are shown 
in figure [1] and given in Table IIVI for Kj^ neu, 
Kl TT^v, Ks irei', K^ irei', K^ tt^v, 
and for the combination. The average, |Kis|/+(0) = 
0.21664(48), has an uncertainty of about of 0.2%. The 
results from the five modes are in good agreement, 
the fit probability is 58%. In particular, comparing 



Heavy Quarks and Leptons, Melbourne, 2008 



3 



0.215 



0.2175 



lavi 

'net 



Kaon WG 



f+(0)xVus 



K 



Le3 



K 



K 



Se3 



K 



K 



e3 



ft3 



0.215 0.2175 
FIG. 1: Display of \Vus\f+iO) for all channels. 

the values of 1/4.(0) obtained from and 
we obtain a value of the SU(2) breaking correction 
'^s;7(2)„p = 2.9(4)% in agreement with the CHPT 
calculation Sgij^2) — 2.36(22)%. Moreover, recent an- 
alyzes on the so-called violations of Dashen's theorem 
in the kaon electromagnetic mass splitting point to 
^su(2) values of about 3%. 

The test of Lepton Flavor Universality (LFU) between 
Ke3 and if ^3 modes constraints a possible anomalous 
lepton-flavor dependence in the leading weak vector 
current. It can therefore be compared to similar tests 
in T decays, but is different from the LFU tests in 
the helicity-suppressed modes 7r/2 and Ki2. The re- 



sults on the parameter r^e 



EjExp 



in 



■SM 



IS 

1.0043 ± 0.0052, in excellent agreement with 
lepton universality. With a precision of 0.5% the 
test in decays has now reached the sensitivity of 
other determinations: ri^eir) = 1-0005 ± 0.0041 and 
133 [2] 



(tt) = 1.0042 ±0.0033 



An independent determination of Vus is obtained from 
Kg2 decays. The most important mode is — > M^^i 
which has been recently updated by KLOE reaching 
a relative uncertainty of about 0.3%. Hadronic uncer- 
tainties are minimized considering the ratio T{K^ — 
^+zy)/r(7r+ -J> ^+z/): 



V„, 



Vud 



K 



X(l +(5cm) 



Using the world average values of BR(_ft'^ A*^*^) 
and of given in Section |TT] and the value of 
r(7r± fi'^iy) = 38.408(7) ^s'^ from we obtain: 
IK. |/|Kd|//f//. - 0.2760 ± 0.0006. 



B. Theoretical estimates of /+(0) and /a /Zit 

The main obstacle in transforming these highly precise 
determinations of |Ks|/+(0) and |Ks|/|Kd|/x//77 
into a determination of \Vus\ at the precision of 0.1%, 
are the theoretical uncertainties on the hadronic pa- 
rameters /+(0) and JkI Itx- By construction, /+(0) 
is defined in the absence of isospin-breaking effects of 
both electromagnetic and quark-mass origin. More 
explicitly /+ (0) is defined by the — > 7r+ matrix el- 
ement of the vector current in the limit = rud 
and Qfcm ^ 0, keeping kaon and pion masses to 
their physical values. This hadronic quantity cannot 
be computed in perturbative QCD, but it is highly 
constrained by SU{3) and chiral symmetry. In the 
chiral limit and, more generally, in the SU{3) limit 
(m,i = md = rUs) the conservation of the vector cur- 
rent implies /+(0)=1. Expanding around the chiral 
limit in powers of light quark masses we can write 
/+(0) — 1 + /2 + /4 + ■ • • where /2 and fi are the 
NLO and NNLO corrections in ChPT. The Ademollo- 
Gatto theorem implies that (/+(0) — 1) is at least of 
second order in the breaking of S'C/(3) This in turn 
implies that is free from the uncertainties of the 
O^p^) counter-terms in ChPT, and it can be computed 
with high accuracy: f-i = —0.023. The difhculties in 
estimating /+(0) begin with /4 or at 0{p^) in the 
chiral expansion. Several analytical approaches to de- 
termine fi have been attempted over the years, essen- 
tially confirming the original estimate by Leutwyler 
and Roos. The benefit of these new results, obtained 
using more sophisticated techniques, lies in the fact 
that a better control over the systematic uncertain- 
ties of the calculation has been obtained. However, 
the size of the error is still around or above 1%, which 
is not comparable to the 0.2% accuracy which has 
been reached for |14s|/-i-(0). 

Recent progress in lattice QCD gives us more opti- 
mism in the reduction of the error on /+(0) below the 
1% level. Most of the currently available lattice QCD 
results have been obtained with relatively heavy pions 
and the chiral extrapolation represents the dominant 
source of uncertainty. There is a general trend of lat- 
tice QCD results to be slightly lower than analytical 
approaches. An important step in the reduction of 
the error associated to the chiral extrapolation has 
been recently made by the UKQCD-RBC collabora- 
tion. Their preliminary result /-i-(0) — 0.964(5) is 
obtained from the unqucnched study with Np — 2 + 1 
flavors, with an action that has good chiral properties 
on the lattice even at finite lattice spacing (domain- 
wall quarks). They also reached pions masses (> 330 
MeV) much lighter than that used in previous studies 
of /+(0). The overall error is estimated to be 0.5%, 
which is very encouraging. 

In contrast to the semileptonic vector form factor, the 
pseudoscalar decay constants are not protected by the 
AdemoUo-Gatto theorem and receive corrections lin- 
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FIG. 2: Results of fits to \Vud\, \Vu^\, and \Vus\/\Vnd\- 



(without imposing unitarity) we obtain; Gckm = 

[IKdl' + IK.I' + iKbl']'^' = (1.1662 ±0.0004) x 
10~^ GeV~^, in perfect agreement with the value ob- 
tained from the measurement of the muon lifetime: 
Gf, = (1.166371 ± 0.000007) x 10"^ GeV"^ The cur- 
rent accuracy of the lepton-quark universality sets 
important constraints on model building beyond the 
SM. For example, the presence of a Z' would affect 
the relation between Gckm and G^. In case of a 
Z' from S'O(IO) grand unification theories we obtain 
mz' > 700 GeV at 95% CL, to be compared with the 
mz' > 720 GeV bound set through the direct collider 
searches Q ■ In a similar way, the unitarity constraint 
also provides useful bounds in various supersymmetry- 
breaking scenarios. 



ear in the quark masses. Expanding //tt in power of 
quark masses, in analogy to /+(0), /k/ fn — 1-1-^2-1-. . . 
one finds that the 0{p'^) contribution r2 is already 
affected by local contributions and cannot be unam- 
biguously predicted in ChPT. As a result, in the de- 
termination of /k/Ztt lattice QCD has essentially no 
competition from purely analytical approaches. The 
present overall accuracy is about 1%. The novelty are 
the new lattice results with Np — 2 + 1 dynamical 
quarks and pions as light as 280 MeV, obtained by 
using the so-called staggered quarks. These analyzes 
cover a broad range of lattice spacings (i.e. a=0.06 
and 0.15 fm) and are performed on sufficiently large 
physical volumes {niT^L > 5.0). It should be stressed, 
however, that the sensitivity of fK/f-rr to lighter pions 
is larger than in the computation of /+(0) and that 
chiral extrapolations are far more demanding in this 
case. In the following analysis we will use as reference 
value the MILC-HPQCD resuff /k/U = 1.189(7). 



C. Test of CKM unitarity 

To determine \Vus\ and \Vud\ we use the value 
\Vus\f+{0) = 0.2166(5), the result iK.I/lKdl/x//. - 
0.2760(6), /+(0) = 0.964(5), and /k/L = 1.189(7). 
From the above we find: |Ks| = 0.2246 ± 0.0012 
from Kes only, and |Ks|/|Kd| = 0.2321 ± 0.0015 
from K£2 only. These determinations can be 
used in a fit together with the the recent evalua- 
tion of Vud from 0+ 0+ nuclear beta decays: 
|y„^|=0.97418±0.00026. This global fit gives Kd = 
0.97417(26) and = 0.2253(9), with xVndf = 
0.65/1 (42%). This result does not make use of 
CKM unitarity. If the unitarity constraint is in- 
cluded, the fit gives = 0.2255(7) and xVndf = 
0.80/2 (67%). Both resuffs are illustrated in Fig- 
ure [21 The test of CKM unitarity can be also in- 
terpreted as a test of universality of the lepton and 
quark gauge couplings. Using the results of the fit 



D. Ki2 sensitivity to new physics 

A particularly interesting test is the comparison of 
the I Vus I value extracted from the helicity-suppressed 
decays with respect to the value extracted from 
the helicity-allowed Kis modes. To reduce theoretical 
uncertainties from fx and electromagnetic corrections 
in K£2, we exploit the ratio BR{Ki2) / BR{-ki2) and we 
study the quantity 



Rl23 



Vus{Kn) Kd(0+ -> 0+) 



Within the SM, i?/23 = 1, while deviation from 1 can 
be induced by non-vanishing scalar- or right-handed 
currents. Notice that in i?j23 the hadronic uncertain- 
ties enter through (/if //7r)//+ (0) . In the case of effect 
of scalar currents due to a charged Higgs, the unitar- 
ity relation between \Vud\ extracted from 0+ — > 0+ 
nuclear beta decays and \Vus\ extracted from Kg'i re- 
mains valid as soon as form factors are experimen- 
tally determined. This constrain together with the 
experimental information of log C'^^'S'S'Af ^^^^ |-,g ^^gg^j 
in the global fit to improve the accuracy of the deter- 
mination of i?;23i which in this scenario turns to be 
^i23lscTiar = 1-004 ± 0.007.^ Hcrc (/k//.)//+(0) has 
been fixed from lattice. This ratio is the key quantity 
to be improved in order to reduce present uncertainty 
on i?;23 . This measurement of i?/23 can be used to set 
bounds on the charged Higgs mass and tan /3. Figure 
[3] shows the excluded region at 95% CL in the Mh- 
tan/3 plane. The measurement of BR(i3 — > tv) can 
be also used to set a similar bound in the M/z-tan /3 
plane. While B ^ tv can exclude quite an extensive 
region of this plane, there is an uncovered region in 
the exclusion corresponding to a destructive interfer- 
ence between the charged-Higgs and the SM ampli- 
tude. This region is fully covered by the K 
result. 
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95% CL from K->nv/ji->nv 
□ 95% CL from B^xv 



200 300 400 500 

charged Higgs mass (GeV/c^) 



FIG. 3: Excluded region in the charged Higgs mass-tan P 
plane. The region excluded hy B ^ tv is also indicated. 



E. A test of lattice calculation 



The vector and scalar form factors f+fi{t) are an- 
alytic functions in the complex t-plane, except for a 
cut along the positive real axis, starting at the first 
physical threshold tth = {"mK + Wtt)^, where they 
develop discontinuities. They are real for t < tth- 
Cauchy's theorem implies that /+.o(0 '^^^ be written 
as a dispersive integral along the physical cut where all 
possible on-shell intermediate states contribute to its 
imaginary part. A number of subtractions is needed 
to make the integral convergent. Particularly ap- 
pealing is an improved dispersion relation recently 
proposed where two subtractions are performed at 
t = (where by definition, /o(0) = 1) and at the 
so-called Callan-Treiman point tcT = {'^'k — wj). 
Since the Callan-Treiman relation fixes the value of 
scalar form factor at tcT to the ratio {Jk / f-n) / f+{^), 
the dispersive parametrization for the scalar form fac- 
tor allows to transform the available measurements 
of the scalar form factor into a precise information 
on (/K//7r)//+(0), completely independent of the lat- 



tice estimates. Figure |4] shows the values for /+(0) 
determined from the scalar form factor slope mea- 
surements obtained using a dispersive parametriza- 
tion and the Callan-Treiman relation, and Jk/ f-ir = 
1.189(7). from result on the FF slope using the disper- 
sive parameterization The value of /+(0) = 0.964(5) 
from UKQCD/RBC is also shown. The NA48 resuh is 
difficult to accommodate. Here one can see that this 
result is also not consistent with the theoretical esti- 
mates of /-i-(0). In particular, it violates the Fubini- 
Furlan bound /+(0) < 1. For this reason, the NA48 
result will be excluded when using the Callan-Treiman 
constraint. 

UKQCD/RBC Pl^vi 



ISTRA+ 



Kaon WG 



W) 



0.9 



0.95 



1 



1.05 



FIG. 4: Values for /+(0) determined from the scalar 
form factor slope using the Callan-Treiman relation and 
fK/h = 1.189(7). 
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